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1  INTRODUCTION 

t 

Over  ihe  past  decade,  finite  element  (FE)  analysis  ol  a  conventionally  launched  projectile  subjected 
to  pressure  loading  has  Income  a  neccssaiy  step  in  the  munition  design  process.  A  methodology  has  been 
developed  which  allows  for  a  two-dimensional,  quasi-static  analysis,  that  is.  a  static  balance  between  peak 
pres.sure  and  a  resisting  body  force  acceleration.  Techniques  for  simplified  modeling  ol  threaded  ^ 

interfaces,  fins,  and  nosetips  have  also  been  formulated. 

Electromagnetic  (EM)  railguns  differ  from  conventional  propellant  gun  systems  m  that  they  do  noi 
use  a  pressure  load  on  the  surface  of  a  projectile  as  the  propulsion  force,  but  a  body  lorcc  acceleration,  t 

provided  by  current  conducted  through  a  solid  armature  in  a  magnetic  field  The  mode  ol  operaii(>n  ol 
EM  guns  has  been  detailed  el.sewhcrc,  and  the  interested  reader  is  encouraged  to  review  ihe  inlormaiion 
presented  in  the  literature  (Price  et  al.  19X7;  Jamison  and  Burden  19X9.  Mongeau  l^»‘Xh 

> 

While  EM  guns  and  gas  pressure  guns  are  different,  the  means  of  examining  a  projeciile N  sirui-iurjl 
integrity  are  the  same.  Various  assumptions  are  required  when  modeling  the  projccnle  geomciry  to  allo« 
for  a  two-dimensional  quasi-static  analysis.  Many  of  the  techniques  vahdaied  in  previous  proicsiilc 
analyses  of  conventional  gun  systems  have  application  in  E.M  projectile  analysis  and  will  be-  uiili/cd  here 
The  details  of  these  modeling  assumptions  and  the  results  of  the  subsequent  analysis  svill  be  a  pnnurv 
focus  of  this  report. 

» 

One  concern  unique  to  the  E.M  projectile  designer  is  the  shonened  rise  time  ol  ihe  lorsc  loading  m 
comparison  to  that  of  gas  pressure.  While  this  capability  may  be  beneficial  by  allow  mg  lor  stioricr  harrei 
lengths,  it  may  possibly  introduce  stress  waves  to  the  structure  leading  to  elevated  stress  states  that  mav 
be  catastrophic,  A  quasi-static  analysis  would  fail  to  account  lor  these  dynamic  cllects  Ihus.  wlien  * 

dealing  with  reduced  rise  times  it  is  imponant  to  perform  a  transient  dynamic  anal \  sis  I  he  commercially 
available  FE  code  .ANSYS  (DcSalvo  and  Ciorman  19X9)  was  used  lor  the  analyses  detailed  in  this  repon 
because  of  its  ability  to  determine  structural  response  under  time-dependent  loading  l  ive  dillereiu  load 

I 

profiles  were  considered  for  the  transient  analysis  The  initial  static  analysis  utili/ed  the  two  peak,  loadings 
from  the  five  load  profiles.  Ultimately,  the  inllucncc  of  pulse  shape  is  assessed  as  a  lunction  ol  the 
projectile  stress-state. 


I 
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It  is  important  to  point  out  that  this  work  was  performed  with  the  intention  of  developing  a  generic 
EM  projectile  design  and  analysis  methodology  using  finite  element  techniques.  The  projectile  modeled 
served  only  as  a  vehicle  on  which  to  refine  the  modeling  technique,  with  the  analysis  not  intended  to 
produce  a  design  with  tactical  applications. 

2.  PROJECTILE  DESIGN  AND  MODELING 

The  projectile  design  under  consideration  is  the  re.sult  of  a  design  study  to  support  the  Cannon  Caliber 
Electromagnetic  Launcher  (CCEML)  effort.  The  methodology  used  to  attain  an  integrated  launch  package 
design  considers  the  target  performance  at  range  and  was  first  used  for  the  0.60-cal.  armor-piercing 
fin-stabilized  discarding  armature  (APFSDA)  projectile  package  (Zielinski  1991).  Subsequent  to  that 
analysis,  the  methodology  was  applied  to  a  cannon-caliber  type  launcher  (Zielinski  1990).  The  design 
procedure  uses  static,  simplified  mechanical  analysis  to  detemrine  rod  and  armature  structural  integrity. 

The  selected  .APFSDA  configuration  which  satisfied  the  CCEML  perlbrmance  requirements  has  a 
2.Vmm  bore  diameter,  a  0.194-lb  (88-g)  launch  mass,  and  a  muzzle  velocity  of  7.710  ft/s  (2,350  m/s). 
Also,  the  analysis  assumed  a  sinu.soidal  current  pulse,  with  the  half-eycle  terminated  at  the  end  of  the  2-m 
barrel  length  (Zielinski  1993). 

A  plastic-composite  fin  as.sembly  with  four  blades  attached  was  incorporated  in  the  design  to  provide 
a  means  of  flight  stabilization.  Figure  1  shows  a  complete  a.ssembly  of  the  EM  projectile.  A  dimensioned 
drawing  of  the  sabot/armaturc  is  shown  in  Figure  2.  Note,  the  integrated  launch  package  is  not 
axisymmctric.  To  simplify  the  analysis  and  keep  it  two-dimensional,  it  was  decided  to  adjust  the  density 
of  the  asymmetric  armature  over  its  rear  0.9  in  (2.3  cm).  The  resultant  density  was  6591  that  of  aluminum 
and  allowed  for  this  rear  portion  of  the  armature  to  be  modeled  as  a  cylinder.  The  specific  calculations 
for  determining  this  value  arc  prc.scnted  in  the  Appendix.  Figures  3  and  4  arc  provided  to  show  the 
geometries  of  the  pcnctrator  and  fin  hub,  respectively. 


A  few  other  assumptions  were  made  to  simplify  the  geometry  of  the  FE  model.  The  mass  of  the  fin 
blades  was  lumped  into  the  fin  hub  to  allow  for  cylindrical  modeling  of  the  tail  section.  A  similar 
approach  was  taken  tn  squaring  off  the  pcnctrator  nose.  In  both  cases,  the  length  of  the  cylindrical  masses 
was  cho.scn  to  keep  the  center  of  gravity  of  the  entire  structure  from  shifting. 
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Initially,  it  was  decided  to  model  each  of  the  34  teeth  of  the  pcnctrator  and  sabot  along  their  interface. 

Figure  5  shows  the  grid  pattern  used  for  each  groove,  with  16  elements  utilized  in  both  penetrator  and 

sabot.  Nodes  were  shared  between  the  pcnctrator  and  sabot  along  the  driving  flank  of  each  tooth  and  are  • 

therefore,  in  effect,  bonded.  A  0.001 -in  (0.00254-cm)  gap  was  modeled  between  the  nondriving  surfaces 

of  the  teeth. 

The  rear  nongrooved  interface  of  the  penetrator  and  sabot  was  modeled  with  gap  elements  provided  * 

in  the  ANSYS  code  which  allow  for  sliding  and  opening  along  the  interface.  It  is  important  not  to  bond 
the  nodes  between  these  two  structures  because  this  would  provide  an  additional  driving  surface  and  allow 
for  stress  relief  of  the  armature  via  additional  load  being  carried  by  the  tungsten  core.  ^ 

The  final  decision  to  be  made  for  the  modeling  effon  concerned  the  way  to  represent  the  EM  loading. 

Previous  analytical  work  has  shown  that  the  current  density  through  the  armanire  is  concentrated  heavily 

towards  the  rear  of  the  armature  (Powell,  Walben,  and  Zielinski  1993).  While  ANSYS  has  the  capability  I 

of  placing  internal  force  loads  within  a  structure,  it  was  decided  to  adopt  a  simplified  approach. 

Therefore,  the  approach  adopted  applied  a  pressure  along  the  back  surface  of  the  armature  to  provide  a 

force  equivalent  to  the  actual  EM  load.  It  was  desired  to  investigate  two  different  load  cases,  one 

providing  a  peak  force  of  39,500  Ibf  (176  kN)  and  a  second  having  a  peak  of  52,600  Ibf  (234  kN).  • 

Newton’s  Second  Law  of  Motion  was  used  to  calculate  accelerations  of  201,000  and  268,000  g’s  for  the 

two  cases.  The  applied  pressures  were  calculated  by  dividing  the  previous  peak  forces  by  the  area  of  the 

armature  over  which  the  pressure  is  to  be  applied.  The  resultant  pressure  values  used  in  the  FE  analysis 

were  70.860  psi  (489  MPa)  for  the  less  severe  load  case  and  94,212  psi  (650  MPa)  for  the  high  load  case.  * 

The  final  modeled  geometry  of  the  projectile  is  displayed  in  Figure  6.  There  are  mote  than  3,300 
elements  due  to  the  detailed  modeling  of  the  grooved  interface  region.  Radial  constraints  were  placed  ^ 

along  the  top  of  the  sabot/annatunc  to  represent  the  presence  of  the  gun  barrel.  This  geometry  model,  with 
the  less  severe  load  case  (176  kN).  was  used  to  begin  the  analysis. 
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3.  RESULTS  OF  THE  PRELIMINARY  ANALYSIS 


The  initial  analysis  of  the  projectile  described  in  the  previous  section  showed  extremely  high  stresses 
in  the  rear  teeth  of  the  penetrator/sabot  interface.  The  model  also  experienced  numerical  convergence 
problems  associated  with  the  use  of  gap  elements.  The  material  properties  of  the  aluminum  were  extended 
to  include  the  plastic  regime  to  allow  deformation  in  an  attempt  to  reduce  stresses.  Also,  the  ANSYS 
manual  suggested  that  ramping  the  load  to  its  peak  value  would  aid  in  the  convergence  of  the  gap 
elements,  and  this  advice  was  followed  (DeSalvo  and  Gorman  1989).  The  mass  of  the  model  was 
calculated  by  the  code  to  be  0.196  lb  (89  g),  which  agreed  with  the  previously  calculated  value  and 
provided  confidence  that  the  geometry  model  of  the  projectile  was  correct  (Zielinski  1990). 

With  multiple  load  steps  and  a  plasticity  model  for  the  aluminum  incorporated,  the  model  was  rerun. 
Again,  the  results  showed  excessive  stress  levels  in  the  rear  teeth,  but  the  sabot  exhibited  severe 
deformation  with  convergence  of  the  gap  elements  achieved.  This  infonnation  led  to  effons  to  reduce  the 
stresses  present  in  the  projectile. 

At  this  point  in  the  analysis  effort,  a  hardware  problem  was  encountered  with  the  computer  system 
in  use.  The  FE  work  had  been  performed  on  a  RISC-based  processor  in  a  HP730  workstation.  This 
system  has  the  ability  to  process  76  million  instructions  per  second  (MIPS).  An  analysis  of  Ihc  EM 
projectile  with  plasticity  material  models,  gap  elements,  and  multiple  load  steps  required  approximately 
8  1/2  hours  to  complete.  With  this  system  down,  the  available  alternative  was  a  Motorola-based  processor 
on  an  Apollo  ring,  capable  of  processing  7.5  MIPS.  Thus,  continuing  on  with  the  present  model  would 
have  taken  at  least  10  times  longer  to  perform  an  analysis  run  than  on  the  HP  machine.  Facing  the 
prospect  of  a  3-1/2  day  run  per  quasi -static  run  on  the  slower  system  and  with  even  lengthier  runs 
resulting  from  transient  analyses,  it  was  decided  to  modify  the  grid  of  the  FE  model. 

Fii.st,  the  grooved  interface  was  modeled  using  a  material  that  represented  the  tungsten-aluminum 
mated  teeth.  This  material  was  given  a  density  equal  to  one-half  the  sum  of  the  den.sities  of  tungsten  and 
aluminum  to  represent  this  interface  to  maintain  a  correct  projectile  mass.  The  material  was  then  given 
an  elastic  modulus  equal  to  that  of  aluminum  since  this  is  the  less  .stiff  of  the  two  materials  which 
comprise  the  interface  region.  This  procedure  has  become  commonplace  in  the  analysis  of  conventional, 
gas  pressure  loaded,  projectile  analysis  (Sorensen  1991).  Adoption  of  these  assumptions  for  the  grooved 
interface  allowed  for  the  use  of  a  much  coarser  mesh. 
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Another  decision  made  to  reduce  the  complexity  of  the  analysis,  and  subsequently  reduce  run  time, 
was  to  cease  with  the  use  of  gap  elements  along  the  nongrooved  penctrator/armature  interface,  thereby 
eliminating  the  need  to  perform  multiple  load  steps.  Instead,  the  elements  in  this  region  were  given  very 
small  values  for  the  clastic  modulus  in  the  radial  direction.  This  was  done  so  as  not  to  allow  for  the 
transfer  of  load  from  the  armature  to  the  penetrator  and  artificially  reduce  the  stress  level. 

The  modifications  previously  outlined  were  incorporated,  the  analysis  was  continued  on  the  Apollo 
ring,  and  the  results  mirrored  those  obtained  previously.  Again,  the  rear  tooth  region  showed  excessive 
stress  levels  with  plastic  deformation  of  the  sabot  leading  to  the  belief  that  the  new  model  on  the  new 
machine  was  consistent  with  the  work  done  previously.  Therefore,  efforts  were  initiated  to  reduce  the 
stress  state. 

4.  ANALYSIS  OF  THE  MODIFIED  PROJECTILE 

It  was  felt  that  increasing  the  length  of  the  grooved  interface  would  greatly  reduce  the  magnitude  of 
the  stress  contours  through  the  projectile.  Therefore,  the  length  of  the  grooves  was  extended  to  span  the 
entire  length  of  the  penctraior/sabot  surface.  The  convention  of  employing  an  adjusted  material  for  the 
annular  grooves  was  continued  and  eliminated  the  need  to  reduce  the  modulus  of  material  along  the  rear 
of  the  interlace  as  had  been  previously  assumed.  The  plasticity  material  model  was  also  dropped. 

The  results  of  this  analysis  showed  that  the  stress  contours  had  been  reduced  to  acceptable  levels.  A 
color  plot  is  shown  in  Figure  7  with  the  stress  levels  listed  at  the  right  and  given  in  pounds  per  square 
inch  (psi).  The  maximum  stress  in  the  rod  is  seen  to  be  L*i7  ksi  (l.OHO  MPa),  well  below  its  yield  value 
of  246  ksi  (1.690  MPa)  A  look  at  the  aluminum  sabot  shows  a  maximum  through  effective  von  Miscs 
stress  of  54  ksi  (.^70  MPa),  also  well  below  iis  yield  value  of  S2  ksi  (565  MPa).  The  stresses  in  the  axial 
direction  were  seen,  as  expected,  in  be  approximately  70  ksi  (4H3  MPa),  the  applied  pressure  load.  The 
maximum  effective  stress  value  does  ni't  have  as  high  a  level  as  the  applied  pressure  due  to  the  radial 
constraints  imposed  on  the  borcriders  which  assists  in  spreading  the  load  to  the  hoop  direction.  An  clastic- 
plastic  analysis  of  the  redesign  was  also  pcriormed  with  little  change  to  the  stress  values,  which  is  as 
expected  since  the  stress  levels  fell  well  short  of  the  malenal  elastic  limits 


Another  ease  which  increased  the  length  ol  the  grooved  interface  over  that  of  the  original  design  but 
left  0.1.3  in  (0.33  cm)  ungrooved  at  the  rear  was  run.  For  this  case,  the  maximum  effective  stress  in  the 
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l-iguiv  7  Slross  ccmunirs  in  Z.Vnun  proicctilc  with,  peak.  U^rcc  (>l  176  kN 

[vnoiraii'r  remained  at  157  ksi  ( l.OSO  Mf’a).  but  the  sabot  through  stress  rose  2()'l  'Urns,  it  was  decided 
let  sla>  'v'.itli  the  I'ulK  grooved  interlace  design  lor  die  high  load  case. 

riie  mote  severe  loail  [irolile  resulted  In  a  [K'ak  stress  ot  210  ksi  (1,44H  MPal  across  the  rod  with  a 
masmium  through  stress  contour  ol  72  ksi  i4db  MPa)  in  the  salxil  A  contour  plot  ol  die  projectile  stress 
levels  is  shown  m  Mgure  K,  with  F  igure  0  providing  a  blown-up  view  ol  the  sabol/amiature.  Again,  these 
are  eltevlive  von  Mises  stress  values  l•..';amination  ol  axial  stress  levels  tilong  tlie  rear  ol  the  annalure 
exceeded  ds  ksi  (he's  \lP;i)  Once  again,  the  explanation  lor  the  difference  between  effective  and  axial 
stress  values  lie^  in  the  boundarv  condition  assumption  of  no  radial  displacement  on  the  bore  riding 
surface  which  assists  in  loading  the  projectile  in  tfie  hoop  direction.  There  is,  how'cver,  a  concern  about 
the  boundarv  coiulition  assumption  employed,  which  is  die  same  as  that  typically  used  in  the  antilysis  of 
gas  pressure  launches  bor  an  17M  launch,  the  magnetic  field  in  the  barrel  tends  to  cause  the  rails  to 
scjiarate  (Zielinski,  Ke/erian,  and  Ik-no  IdSS),  thus  freeing  the  annature  to  move  radially  outward 
Another  point  to  note  trom  f-igure  d  is  the  large  compression  stress  along  the  base  of  the  armature  at  the 
rear  ol  the  interlace  Pan  of  the  reason  for  this  localized  high  stress  is  the  s(|uare  geometry  ot  the  model 
III  lliis  region  as  well  iis  the  verv  small  element  size  used  to  model  the  geometry  This  area  ot  stress 
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concentration  is  something  which  could  be  handled  by  a  change  in  geometry  of  the  design.  Since  the 
scope  of  this  effort  is  on  development  of  the  analysis  techniques  and  not  on  design  techniques,  no  physical 
changes  were  made.  In  an  actual  design  effort,  this  stress  concentration  would  need  to  be  further 
addressed. 

Another  case  of  the  high  load  profile  was  run  with  the  radial  boundary  constraint  removed,  providing 
the  armature  with  complete  freedom  of  outward  movement.  While  it  is  recognized  that  this  too  is  not  an 
accurate  physical  model  (the  rails  do  provide  some  restriction  on  the  armature’s  radial  motion),  it  was  felt 
that  this  case,  accompanied  with  the  previous  one,  would  bracket  the  range  of  possible  armature  motion 
and  the  stresses  which  result.  Lacking  a  true  representation  of  the  rail  separation,  this  was  felt  to  be  the 
best  approach. 

The  results  of  the  unconstrained  case  are  shown  in  Figure  10,  and  not  unexpectedly,  the  effective 
stress  through  the  armature  has  risen  appreciably  to  92.5  ksi  (638  MPa),  a  27%  increase.  The  lack  of  a 
radial  constraint  reduces  the  load  carrying  capability  in  the  hoop  direction.  This  also  causes  an  increase 
in  the  penetrator  stress  to  223  ksi  (1,538  MPa).  This  then  points  up  the  importance  of  establishing 
realistic  boundary  conditions  when  performing  the  FE  analysis.  A  poor  choice  of  boundary  conditions 
can  lead  to  a  projectile  which  lacks  sufficient  robusmess  to  withstand  launch  conditions  being  deemed 
acceptable.  This  highlights  the  point  that  the  FE  codes  are  tools  which  provide  results  only  as  accurate 
as  the  physical  reality  of  the  input. 

5.  TRANSIENT  ANALYSIS 


Previous  work  has  shown  the  shorter  rise  times  typical  of  EM  launches  can  lead  to  an  amplification 
of  the  stress  levels  due  to  the  generation  of  stress  waves  (Bannister,  Burton,  and  Drysdale  1991).  While 
that  work  focused  on  a  lengthier  projectile  subjected  to  a  much  shorter  rise  time  loading  than  the  cases 
under  consideration  here  for  the  23-mm  projectile,  it  remains  important  to  consider  the  effects  of  stress 
waves  for  projectile  loadings  whose  rise  times  approach  the  time  required  for  these  waves  to  sufficiently 
dampen  out  A  saboted  projectile  is  an  undamped  system,  thus  requiring  a  stress  wave  to  traverse  the 
length  of  the  penetrator  numerous  times  before  attaining  a  steady-state  value  approaching  zero.  Work 
done  investigating  stress  waves  uirough  bare  rods  showed  that  rise  times  greater  than  five  times  the  transit 
time,  the  time  required  for  a  pulse  to  travel  down  the  length  of  the  penetrator  and  back,  were  of  sufficient 
duration  to  ensure  sufficient  damping  of  the  stress  waves  (Bcrxlcr  1989). 
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Table  1.  Acceleration  Profile  Parameters 


Acceleration 
Case  No. 

Peak  Acceleration 

(g) 

Time  of  Peak 
Acceleration 
(ms) 

Time  of  Travel  to 
Muzzle  exit 
(ms) 

1 

268.000 

0.80 

3.55 

2 

268,000 

1.45 

8.025 

3 

201,000 

0.45 

2.70 

4 

201.000 

0.80 

5.50 

5 

201.000 

1.45 

4.25 

The  acceleration  profiles  were  resolved  to  pressure  values  every  25  ps,  which  deteimined  the  number 
of  ramped  load  steps  used  in  the  FE  run.  For  instance,  case  3  had  a  muzzle  exit  time  of  0.0027  s,  which 
results  in  108  load  steps  (0.0027/0.000025).  The  other  important  consideration  for  the  transient  analysis 
was  selection  of  an  integration  time  step  which  is  established  by  the  selection  of  the  number  of  iterations 
per  load  step.  The  ANSYS  code  suggests  guidelines  to  ensure  the  integration  time  step  is  sufficiently 
small  to  resolve  the  motion  of  the  structure,  as  well  as  account  for  possible  wave  propagation  (Desalvo 


and  Gorman  1989)  L  sing  the  ANSYS  criteria,  a  value  of  0.1  ps  was  chosen  as  the  integration  time  step. 
This  then  provided  lor  250  iterations  per  load  step  in  the  FE  code. 


Having  deemed  a  vcr\  small  time  step  necessar>'.  it  was  decided  to  derive  another  FE  model  of  the 
projectile  geomeirv  uiili/ing  a  coarser  grid  to  reduce  the  time  required  to  perform  a  single  analysis  run. 
The  model  was  reduced  from  7.5.5  elements  to  287  elements  and  rerun  with  a  quasi-static  pressure  load. 
There  was  no  appreciable  difference  in  the  stress  contours  between  these  two  models,  so  the  more  coarsely 
gridded  projectile  was  earned  forth  for  use  in  the  transient  analysis. 

The  five  transient  cases  were  run  with  a  radial  constraint  applied  along  the  entire  boteriding  surface 
of  both  the  sabot  and  armature.  Stress  plots  of  the  projectile  at  the  time  of  maximum  equivalent  stress 
are  shown  in  Figures  12- lb  for  the  five  acceleration  loading  profiles.  Neither  the  sabot  nor  pcneirator 
reach  values  exceeding  their  respective  yield  values  The  results  of  the  transient  analysis  show  agreement 
w  ithin  2fr  of  the  calculated  stress  values  for  the  quasi-static  analysis  runs.  Again,  it  is  noted  that  removal 
of  the  radial  displacement  constraint  would  likely  result  in  stresses  exceeding  the  saboi  matenal's  yield 
strength. 

Depicted  in  Figures  lb-21  arc  plots  of  the  stress  time  history  of  an  element  located  in  the  center  of 
the  penetraior.  These  plots  rellcci  the  various  rise  and  fall  times  of  the  acceleration  profiles  used  as  input. 
It  is  noticed  that  every  trace  exhibits  an  o.scillaiory  respoasc  following  the  attainment  of  peak  stress.  To 
determine  whether  this  is  a  transient  effect  due  to  the  unloading  of  the  current  through  the  projectile  or 
an  artifact  of  clement  hour-glassing,  it  is  necessary  to  calculate  the  transit,  time,  T.  This  transit  time  is 
defined  by 

.p  _  total  pcneirator  length  traversed 
wave  speed  of  the  matcnal 


The  wave  speed  for  a  material  is  expressed  as  y^E/p  ,  where  E  is  the  clastic  modulus  of  the  material 
and  p  is  its  density.  For  the  tungsten  pcneirator,  this  results  in 


T  = 


(2)  (5.865  in) 


_  =  47.5m  s 


(44x10  psi 


/  f  0.658  Ibm/in  ^ 

ISI)  /  I - 

'  586.4 


1618 
25221 
48823 
72426 
96028 
119631 
143233 
166836 
190438 
214041 


STLP=59 
ITER=250 
TIME  =  0.  00145 


[~1 

Ft.5v~iS1_ 


CZD 


1608 

25230 

48852 

72474 

96096 

119718 

143339 

166961 

190583 

214205 


STEP-Xd  .  .  r 
TIHE-0445«Er03 


3d4iO 

54X22 

7X764 

60406. 

X07048 

124690 

142332 

15997,4 


4TEP^33lii« 
ITER-250  " V 
TIME-0.  BOOE-t 

1205/ -- 

n  36562s 
S  54240 
S  71919  ; 
S  69597^- 
=  107.275 
=  124954 
S  142632 
™*  160310 


74 


iltP=59 
I TER-250 
TIME-0.  00145 


1207 

18949 

36691 

54433 

72175 

89917 

107659 

125401 

143143 

160885 


25 


Pressure  (psi) 


Figure  20.  Stress  history  of  pcnetrator  clement  subjected  to  acceleration  case  4. 


Figure  21.  Stress  history  of  pcnetrator  clement  suhiected  to  acceleration  case  5. 
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A  comparison  of  the  penetrator’s  calculated  wave  transit  time  vs.  the  response  frequency  of  the  rod  seen 
in  Figures  16-21,  shows  that  the  oscillations  are  much  slower  than  the  transit  time  and  arc  thus  not 
attributable  to  the  transient  dynamics  of  the  unloading  of  the  current  pulse.  It  is  believed  that  the 
fluctuations  in  stress  levels  are  a  result  of  element  hour-glassing  effects. 

6.  CONCLUSIONS 

This  work  laid  out  a  methodology  for  the  analysis  of  EM  projectiles  using  conventional  FE  techniques. 
A  pressure  loading  on  the  rear  of  the  armature  was  utilized  to  simulate  the  magnitude  of  the  true 
electromagnetic  body  force.  An  improvement  to  the  technique  would  be  to  incorporate  the  FE  capability 
of  applying  forces  directly  to  internal  nodes  of  the  modeled  structure.  There  are  currently  several  codes 
written  or  under  development  which  predict  the  current  distribution  through  the  armature 
(Powell  et  al.  1993)  whose  results  could  be  applied  as  input  loads  for  the  structural  analysis.  However, 
the  technique  as  demonstrated  here  provides  a  relatively  good  approximation  of  the  loading  conditions  and 
is  rather  easy  to  implement. 

The  analysis  effort  points  out  the  importance  of  applying  realistic  boundary  conditions.  Utilizing  a 
radial  boundary  constraint  as  typically  done  for  the  analysis  of  projectiles  in  conventional  gun  systems  can 
lead  to  a  faulty  design  being  deemed  acceptable.  Obviously,  such  an  erroneous  conclusion  could  be 
catastrophic.  The  separation  of  the  rails  inherent  in  the  EM  guns  is  of  concern  to  projectile  designers  and 
should  be  addressed  in  future  analysis  work. 

Finally,  while  the  transient  analysis  performed  on  the  projectile  during  this  investigation  produced  little 
difference  compared  to  the  quasi-static  cases,  it  is  important  that  the  EM  projectile  designer  go  through 
the  process  of  examining  the  dynamic  loading  effects.  This  becomes  especially  true  as  the  rod  length  is 
increased  and  the  rise  time  is  reduced,  resulting  in  insufficient  time  for  the  oscillations  of  any  stress  waves 
to  decay  away.  Longer  rods  mean  increased  transit  times  which  may  lead  to  amplification  of  stresses  as 
they  approach  the  rise  time  of  the  acceleration  profile.  While  Bender’s  work  with  bare  rods  has  shown 
rise  times  in  excess  of  five  times  to  be  of  sufficient  length  to  ignore  the  uansient  wave  effects,  others 
working  with  saboted  projectiles  typically  use  a  factor  of  ten  between  the  rise  time  and  transit  time  before 
ignoring  the  wave  phenomenon  (Hopkins  1993). 
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The  EM  projectile  analyzed  in  this  investigation  does  not  have  an  axisymmetric  cross  section  as  seen 
from  Figure  2.  It  was  desired  to  avoid  having  to  use  a  three-dimensional  element  model  to  simplify  both 
the  construction  of  the  projectile  geometry  and  the  time  required  to  solve  the  problem.  The  adopted 
approach  was  to  model  the  rear  armature  shown  below  (Figure  A-1)  as  a  cylinder  with  an  equivalent  mass. 
Since  only  the  rear  0.9  inch  (2.29  cm)  of  the  armature  is  replaced  by  a  cylinder  of  the  same  length,  there 
is  no  need  to  worry  about  a  shift  in  the  center  of  gravity. 

The  First  step  was  to  calculate  the  weight  of  the  actual  armature  configuration.  This  required  a 
determination  of  the  volume,  and  in  turn,  the  area  of  the  cross  section.  The  area  was  determined  by 
subtracting  the  area  of  the  penetrator,  and  the  area  of  above  the  flat,  from  the  area  of  the  full 
cylinder,  Aj^u,  of  radius  0.452  in  (1.15  cm).  The  areas  are  calculated  as 

.  0.321  in^ 

A™.  ^  ^  .  0.030  in^ 


A«.  =  i  (rad  A  -  sin  A)  -  ^ (0.452  in)^  (1.92  -  0.94)  ^  0.100  in^  . 

*  2  2 

where  rad  A  is  the  radian  measure  of  angle  A  shown  in  Figure  A-1  (Marks  1951).^  Thus, 

^Ktu«i  "  |Afuii  A^  -  A^j]  .  (A-1) 

where  L  is  the  0.9-in  length  and  pj(.,u4|  is  the  density  of  aluminum,  0.98  Ib/in^.  Substitution  of  values  into 
equation  A-1  yields  a  ma.ss  of  0.0168  lb  (0.00076  kg).  Division  of  this  mass  by  the  volume  of  the  full 
cylinder  minus  the  rod  volume  will  provide  a  density  for  the  modified  armature  geometry.  That  is. 


Marks,  L.  S.  Mechanical  Engineers*  Handbook.  New  York:  McGraw-Hill  Book  Company,  Inc..  1951. 


» 


P  fudge 


mass 

(Afuu  A^)  L  ■ 


(A-2) 


> 


Again,  substitution  of  values  into  equation  A-2  gives  a  value  of  0.064  Ib/in^  (1,770  kg/m^)  which  is  used 
for  the  armature  material  in  the  geometry  model. 
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